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Abstract 

In this work we explore 1 + 1 dimensional p-wave superconductors using the probe D-brane construction. 
Specifically, we choose three intersecting D-brane models: D1/D5, D2/D4 and D3/D3 systems. According 
to the dilaton running behavior, we denote the former two systems as nonconformal models and the last 
system as conformal. We find that all three models are qualitatively similar in describing superconduct- 
ing condensate as well as some basic features (such as the gap formation and DC superconductivity) of 
superconducting conductivity. There also exist some differences among the three models as far as the AC 
conductivity is concerned. Specifically, for D3/D3 model there is no peak at nonzero frequency for the 
imaginary part of the conductivity, which is present in the nonconformal models; their asymptotic behaviors 
are different — for D1/D5 the real part of the AC conductivity approaches one at large frequency limit, for 
D2/D4 it slowly goes to a certain nonzero constant smaller than one and for D3/D3 it goes to zero. We find 
that the profile of the AC conductivity for the D1/D5 system is very similar to that of higher dimensional 
p-wave superconductors. 
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I. INTRODUCTION 



The applications of AdS/CFT correspondence JT), more generally gauge/gravity duality, to in- 
vestigations of strongly coupled system have gained broad interest varying from QCD phenomena 
at low energy to strongly correlated condensed matter physics, see, e.g., § for recent reviews. One 
of the most interesting applications is the construction of superconducting-like phase transition. 
Following the pioneering work of ||], holographic superconductors have been constructed in 
by putting the Abelian Higgs model or SU(2) gauge field into the AdS black hole geometry. When 
the Hawking temperature is decreased to some critical value, the black hole background becomes 
unstable against perturbations and gets hair by condensing some field in order to cure the insta- 
bility. This can be considered as holographic realization of the superconducting phase transition. 
This kind of construction of holographic superconductors takes the (asymptotically) AdS black hole 
spacetime as the starting point. In some sense, this construction should be taken as a bottom-up 
approach to the gravity dual of strongly interacting superconductor as the theory in the bulk is 
directly written down from phenomenological point of view. Recently, there appear some works on 
the UV completion of these phenomenological superconducting models, see, e.g., (7j, by embedding 
the holographic superconductors into the superstring/M-theory or gauged supergravity. Another 
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top-down approach to holographic superconductors where the dual field theory is known is based 
on probe D-brane in a black p-brane supergravity geometry. Such a holographic superconductor 
has been established in || where there is a p meson condensate. Moreover, a stringy mechanism 
for the condensation process has been described. For the Sakai-Sugimoto model M, a holographic 



superconductor involving a p meson condensate has been described in [10], based on earlier results 
on p meson condensation in [11]. In [12] the D3/D5 system, constructed to be a defect theory, was 
used to study different condensates corresponding to vector/scalar modes on the flavor D5 brane 
world volume. 

The models mentioned above are all concerned on higher dimensional spacetime, say planar or 
1+3-dimensional superconductors. However, 1+1 dimensional system is also of great importance 
and interest in condensed matter physics. It is therefore of large interest to see what holographic 
methods tell us about such systems. The work of [13] takes the D3/D3 brane system (This model 
was first studied in [O]) to model 1+1 dimensional strongly coupled quantum liquid and found 
some interesting properties different from higher dimensional counterparts. Similar studies taking 
the BTZ black hole geometry as the gravity background can be found in |15|]. Some aspects of 



holographic Luttinger theorem have been discussed in detail in [16|. More recently, holographic 
s-wave superconductors in 1+1 dimensions have been constructed and well-studied in |l7|-p^|l by 
introducing Maxwell- Scalar system into the BTZ black hole geometry. Other recent works for 1+1 



dimensional boundary theory from holography can be found in [20-22]. In |23] local aspects of 



1+1-dimensional superconductivity, which followed the pioneering work of [24|, have been inves- 
tigated in a bottom-up approach. This ensures a true superconductor where a local symmetry is 
spontaneously broken instead of a superfluid where a global symmetry is broken. 

In this paper we take the probe D-brane approach to study some aspects of p-wave superconduc- 
tor in 1+1 dimensional spacetime. More specifically, we take three intersecting D-brane systems as 
our starting point: D1/D5, D2/D4 and D3/D3. For stability of these brane systems, we keep them 
supersymmetric in the sense that the ND number of each system is 4. For simplification, our study 
is limited to the zero quark mass and probe limit (i.e., the backreaction of the flavor probe brane 
on the background geometry is neglect). As in ||, we embed two coincident probe D-branes into 
the black p-brane geometry and just reserve the Yang-Mills truncation of nonlinear DBI action 
for probe D-brane. As mentioned earlier, one advantage of taking the intersecting D-brane system 
is that its field theoretical side is known, here they are supersymmetric gauge theories coupled 
with fundamental matters. The superconducting condensates are nonabelian gauge fields, living 
on the probe flavor branes, and have one Lorentzian index, indicating the phase transition is of the 



3 



p-wave type. Our numerical results for the D3/D3 model are very similar to those of p3| . However, 
one main advantage of the D3/D3 system is that its dual field theory is known (the defect 1+1 
dimensional CFT) and the condensed operator can be explicitly written down, see eq. (p^). 

Although these probe D-brane constructions for p-wave superconductor seem very general as 
well as similar for different D-branes, the results are not always the same. Specifically, we find 
that nonabelian condensates for the three models are qualitatively the same, they have mean 
field behavior near the critical temperature and approach to some fixed constants at very low 
temperature. When turning to the electromagnetic response, we find that, for the real parts of the 
conductivities, the three models give some features in common, like the DC infinite conductivity 
and the gap formation. However, for D3/D3 model, there is no peak at nonzero frequency for 
imaginary part of the AC conductivity compared to the others. Actually, the conductivity formula 
for the D3/D3 system is very different from the other two models due to the nontrivial behavior 
for the electromagnetic fluctuations near the AdS boundary. Another main difference among these 
models is that their asymptotic behaviors of the AC conductivity are quite different. The real part 
of the AC conductivity approaches one at large frequency limit for D1/D5 system and goes to a 
nonzero constant (much smaller than one and also very slowly) for D2/D4 model while for D3/D3 
system it tends to zero. 

The rest of this work is organized as follows. In section |n| we shortly review previous studies on 
quantum field theory in 1+1 dimensional spacetime by taking the approach. We then introduce the 
models studied in this work and numerically solve the nonlinear equations of motion for background 
fields. With these numerical solutions, we plot the superconducting condensate as well as the 
free energy versus the dimensionless temperature. Section |l | is devoted to the study of the 
electromagnetic fluctuation. We plot the AC conductivities for all the three systems and give some 



comments on the results. We end with a short summary and some discussions in section IV 



II. FLAVOR P-WAVE SUPERCONDUCTORS FROM INTERSECTING D-BRANES 
A. Some reviews 

In this subsection, we shortly review previous studies about 1+1 dimensional quantum field 
theory by taking the AdS/CFT approach. The first aspect we intend to mention here is about 
holographic quantum liquids in 1+1 dimensions, initially studied in [13, 15]. Although the authors 



in |T^] take defect D3/D3 intersecting D-brane models as the starting point, their results are in 
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quite agreement with those of [15], which directly takes the charged BTZ black hole to compute 
correlation functions for probe scalar, spinor and vector operators in this geometry. In contrast 
with higher dimensional defects, a persistent dissipationless zero sound mode is found in [|l3|]. The 
correlation between log periodicity and the presence of finite spectral density of gapless modes is 



seen in [15]. Meanwhile, the real part of the conductivity (given by the current-current correlator) 
also vanishes as u — > as expected. The fermionic Green's function shows quasiparticle peaks with 
approximately linear dispersion but the detailed structure is neither Fermi liquid nor Luttinger 
liquid and bears some similarity to a " Fermi-Luttinger" liquid. As will be seen in later sections, to 
some degree, our numerical results for the AC conductivity in D3/D3 model is in agreement with 
these conclusions. 

The second point is about holographic realization of symmetry breaking in 1+1 dimensional 
spacetime, which has been discussed in [|17], 18]. The authors in [18] considered a system of 3D 
gravity coupled to matter to study the symmetry breaking phases in 1+1 dimensional spacetime. 
To be specific, they model symmetry breaking phases of a strongly coupled 1+1 dimensional CFT 
as black holes with scalar hair. It concluded that, in the case of a discrete symmetry, these theories 
admit metastable phases of broken symmetry. Moreover, the 3D Einstein-Maxwell theory shows 
continuous symmetry breaking at low temperature. The latter conclusion can be used to con- 
struct holographic s-wave superconductors in 1+1 dimensional spacetime. Intuitively, continuous 
symmetry breaking in 1+1 dimensional spacetime seems to contradict with the Coleman-Mermin- 
Wagner (CMW) theorem [^5|] which states that in 1+1 and 2+1 space-time dimensions at finite 
temperature, spontaneous continuous symmetry breaking is impossible. However, in the large N 
limit, these lower dimensional systems can have another phase in which the continuous symmetry 
is almost spontaneously broken and the fall of the correlation functions is of the power law type as 
pointed out in |26|]. Under the framework of the AdS/CFT correspondence, one may expect that 
0(1/N) corrections will wash out the symmetry breaking phase. Actually, this indeed happens 



as was explored in [27] for the AdS^ black hole. It is then reasonable to believe that this idea 
also holds in the three dimensional bulk case. Given these facts, it is very interesting to see what 
holographic method can tell us about the superconducting phase transition in two dimensional 
spacetime with vector order parameter. What we found is that, when working in the large N limit, 
which is the fundamental assumption of AdS/CFT correspondence, holographic superconductor 
symmetry breaking in 1+1-dimensional spacetime does happen. 

Another interesting feature for 3-dimensional bulk theory is the chiral anomaly, which has been 
studied in [21] with applications to condensed matter physics and in [22] for correctly producing 



meson spectrum. In the bulk side, the chiral anomaly is the Chern-Simons term in the bulk 
action. The author in [^T| studied the holographic description of finite-density systems in two 
dimensions. Quite interesting, it was shown that the chiral anomaly for symmetry currents in 2 
dimensional CFT completely determines their correlators. The important exception is a CFT with 
a gauge theory to which we may couple an external current, as in the probe D3/D3 system or the 



putative dual to the charged BTZ black hole. In [22], the defect D2/D8 brane model was used 
to holographically realize large Nc massless QCD in two dimensional spacetime. The flavor axial 
anomaly is dual to a three dimensional Chern-Simons term which turns out to be of leading order, 
and it affects the meson spectrum and holographic renormalization in crucial ways. It was also 
shown that an external dynamical photon acquires a mass through the three dimensional Chern- 
Simons term as expected from the Schwinger mechanism. Massless two dimensional QCD at large 
Nc exhibits anti-vector-meson dominance due to the axial anomaly. Explicitly, one still cannot 
give a general argument on the effect of the chiral anomaly term on the dynamics. However, the 
models studied in this work do not contain the chiral anomaly terms, which greatly simplify our 
analysis as well as numerical computations. If we go beyond the massless limit for the flavor quark, 
the Chern-Simons (corresponding to chiral anomaly) will appear and we leave the study along this 
line for future work. 

Before concluding this subsection, we briefly comment the models used in this work. Although 
the three models considered here look quite similar, there are still some differences among them, 
which result in different features for the AC conductivity, which we will reveal in later sections. The 
D1/D5 and D2/D4 models are nonconformal in the sense that the dilaton profiles are non-trivial. 
Then, according to the AdS/CFT correspondence, the gauge coupling constant on the dual field 
theory side run along the radial direction, which can be thought of as the energy scale. We may 
naively denote these two models as the nonconformal SYM gauge theory coupled to fundamental 
matter. However, the D3/D3 model is a conformal one. These differences together with the specific 
background spacetime tell us that the gauge field living on the probe D-brane has quite different 
asymptotic behavior (here we mean the behavior at the AdS conformal boundary) for these models. 
In particular, there is a logarithmic term for the boundary expansion of gauge field in the D3/D3 



case, which will be clear in II C. This makes one to identify the coefficient of the logarithmic term 

instead of the constant term as the source. This never happens for higher dimensional (asymptotic) 

AdS geometry casd]]. Furthermore, due to this point, one have to give a new prescription for the 

1 However, this happens in choosing the Lifshitz black hole geometry to construct the strong coupling superconductor 
with dynamical exponent, as seen in |28|. 
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computation of the transport coefficients in D3/D3 system. The details on these issues will be 
given in later sections. In [||, |To|, the authors mainly focused on the decoupled sector^], which has 
similar equation of motion as in the s-wave case. However, here we directly encounter the coupled 
modes, and need to define a gauge invariant variable to plot the conductivity. This makes our 
numerical computations complex, especially for the D3/D3 system. 



B. Equation of motion for the background fields 

The gravity dual of p-wave superconductor was first constructed in [||, || by putting the SU(2) 
gauge field into the AdS^ black hole geometry. The action for this gravity system is 

S = 2^ / d ' x ^[R - 2A - \f; v F^I (1) 

where the field strength tensor of nonabelian SU(2) gauge field A°^ is defined as F£ u = d^A^ — 
d v At + e abc A b ^A^ with totally antisymmetric tensor e 123 = +1. Working in the probe limit allows 
us to ignore the backreaction of the gauge field on the background geometry. Therefore, we can 
fix the black hole geometry as the Schwarzschild-AdS spacetime and study nonabelian gauge field 
in this curved geometry. The anisotropic features of the p-wave superconductor were explicitly 
reflected on different behaviors of the conductivities along x and y directions. Followed by this 
work are some developments of holographic p-wave superconductors, sec [[3(J for an incomplete list. 

However, all of these works are taking higher dimensional AdS black hole geometry as the 
starting point. We here use the probe brane method, which is first applied to holographic p-wave 
superconductor in ||, to explore some properties of 1+1 dimensional p-wave superconductors. We 
choose the p-brane geometry as our background metric, 

,2 



dsl = H ~ 1/2 (~f(p)dt 2 + dx 2 ) + H 1 ** (jfa + p'dnl^ 



(2) 



with 



e* = H^, H=^\ /(P) = l"(^) 7 " (3) 

where x — (x j X j * * * . X p ) denotes the p-dimensional space of the black p-brane. Explicitly there 
is a horizon in above metric at p = po and the Hawking temperature is 

7- P fPo\^ (4) 



4vrL \L 



The coupled sector in the probe D-brane setup has been revealed recently in 
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For our purpose, we will choose p = 1, 2, 3 and parameterize the internal spaces £l$- p as follows, 

dfi? = d6 2 + cos 2 6dSl + sin 2 6dS' 3 2 , 
dn 2 6 = dO 2 + cos 2 8dS 2 + sin 2 ddS$, 

dn 2 = d6 2 + cos 2 ddi 2 + sin 2 9dS 2 . (5) 

The probe D-brane extends along (t, x , p, S3) for D1/D5 model, (t, x , p, S2) for D2/D4 model 
and (t, x l ,p, for D3/D3 system, respectively. The embedding profile for the probe D-brane 
can be parameterized by 6{p). As mentioned before, we here consider the zero quark mass case, 
i.e., 9{p) = and leave the effect of the nonzero quark mass for future research. With these 
assumptions, the induced metrics on the flavor probes are 

44 = {-f(u)dt 2 + dx 2 ) + ^ + %udS 2 , f(u) = l-u e , 

ds 2 = ) 5/2 u- 5 /2 (-f(u)dt 2 + dx 2 ) + ^n" 3 / 2 ^ + ^uVHSl f(u) = 1 - n 5 , 

Pq n ' Pa 

d4=(f) 2 l?{-f(u)dt 2 + dx 2 ) + £j^ + L 2 de, f{u) = l-u\ (6) 

where in above formulae we have transformed holographic coordinate p to a finite interval [0,1] 
by transformation po/p = u because we found it is more convenient to work with u coordinate 
for numerical computations. In this new coordinate system, u = denotes the AdS conformal 
boundary where the dual field theory lives and the horizon is located at u = 1. 

Two coincident probe D-branes in above black hole geometries have U(2) gauge symmetry on 
its worldvolume. However, we in this work concentrate on its nonabelian subgroup SU(2) for the 
purpose of inducing a p-wave superconducting phase transition. In the probe limit, dynamics of the 
probe Dq-brane is fully determined by nonabelian DBI action^ and we only reserve the Yang- Mills 
truncation of it, 

s = _Mi J d^xe-^gF^^ (7) 

with the determinant g calculated from the induced metric in eq. (|6|) . One further assumption which 
will simplify our computation is that we do not consider the internal coordinates dependence of 
the SU(2) gauge field. Then we can integrate out the internal space in eq. (0) 

S = -N q J d 3 xV^GF* u F aflu , (8) 

3 Note that, for the intersecting D-brane models considered here, there is no Chern-Simons term contribution to the 
D-brane action. For one thing, we take the zero mass limit for the probe D-brane embedding profile; for another, 
we do not consider the excitation of the gauge field along the internal space. 
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where \[—G = ^J—ge ^gs with g$ the internal metric and Af q is a model dependent factor, which 
is irrelevant for later computations. The equation of motion from this action is of the form, 

d^V^GF^] + V^Ge abc A b ^F c ^ = 0. (9) 

We consider the chemical potential induced superconducting like phase transition. To achieve 
this goal, we turn on one time component of nonabelian SU(2) gauge field, say Aq(u) ^ (here, 
we only consider homogenous holographic superconductor, i.e , the background ^4g and A\ only 
depend on holographic coordinate u.). Similar to the arguments in S, one can show that when the 
chemical potential (whose meaning will be clear later), provided by the nonnormalizable mode of 
Aq, is increased to some critical value the black hole will get unstable against perturbations. This 
instability can be cured by condensing some component of the gauge field A 1 ^, say ^4}(n). More 
specifically, we consider following ansatz^|for the hairy black hole, 

A = A%r 3 dt + A\r l dx. (10) 

As emphasized in the first two references of ||, the solution with nontrival profile for A\ will be 
the new ground state when T < T c and this new ground state can be interpreted as a p meson 
superfluid. In what follows we take the D3/D3 system as an example to explicitly write down the 
condensate operator from field theory point of view. More detailed discussions can be found in the 
first two refs of ||. The isospin chemical potential, provided by nonnormalizable mode of Aq, is 
introduced as source of the operator 

J 3 oc Vw 3 7oV' + <pcr 3 do(p (11) 

where ip = (tpui^d) arid 4> = {4>u,<l>d) are fundamental quarks and squarks under gauge group 
SU(2); a % denotes the Pauli matrices and 7^ the Dirac matrices in two dimensional spacetime. In 
the same way, the condensate operator, holographically provided by normalizable mode of bulk 
field A\, takes the following form 

J\ oc VidSiV 1 + ^dxcj). (12) 

The condensate in eq.(|io|) breaks both the SU(2) and translational invariance in the bulk com- 
pletely. However, one should keep in mind that the broken symmetry considered here is the flvaor 
symmetry, like that of the QCD theory, which is somewhat different from the symmetry being bro- 
ken in studies such as ||, [|. The symmetry breaking pattern going through the superconducting 

4 Another configuration for the background has been considered in ]B|, ^| . But it is unstable and have a higher energy 
compared to the one in eq. (lid). 
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phase transition can be understood in the following way. The nonzero value of Aq at the AdS 
conformal boundary explicitly breaks the SU(2) to its subgroup £7(1)3, generated by rotation in 
the colored 12-plane. We can identify this residual unbroken symmetry as the electromagnetic 
symmetry. Strictly speaking, this identification is not right because gauge symmetry in the bulk 
corresponds to a global symmetry on the boundary field theory side. However, this model can 
produce many superconductor-like features. We therefore ignore this difference and just go ahead. 
The U(l) symmetry should be spontaneously broken when going through a superconducting phase 
transition. It is implemented by the nonzero expectation value for the operator O dual to A\. 
Additionally, we need impose the source for O to be zero for spontaneous symmetry breaking. 
Plugging the ansatz fllOD into the equation of motion (|9|) results in 



u 



V-- 

u f(u) 

/'(«) _ 1 
u 



/(«) 



^<j> = 0, 



(13) 



P{u) 







where = (^-O'^-i) f° r D1/D5 model. Similar results for D2/D4 and D3/D3 systems are 

listed as below, 



M = 0, 



T /(«) 



(14) 



/(«) 



<j> 2 i}) = 



with (0, iji) 



AttT 



{A%A\) for D2/D4 and 



1 

+ -< 
u 



/(«) 



f{u) | 1 
f(u) u 



1 



(15) 







with (4>,ip) = rjp (^4q,^4^) for D3/D3 one. In above equations, the prime denotes derivative with 
respect to u and this notation convention will be used in later representations. 



C. Solution for the background fields 

Due to the nonlinear coupling between <fi and ifi in eqs. (|j3|,14jL5|), we turn to a numerical 
shooting method to solve them. The philosophy of our numerical approach is that we first find a 
power series solution for ip and <j> near the horizon. We then take these near horizon solutions as 
initial conditions to numerically integrate these fields from the horizon to the conformal boundary. 
To the conformal boundary (here, represented by u = 0), we impose that the source for the operator 
O is zero. This condition can filter out the wanted solutions. 
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We now have a look at the asymptotic behavior of the background fields <fi and ip. Near the 
conformal boundary, we can get the following asymptotic behaviors from Frobenius analysis of 
eqs. (13,H|Jl5D near the singularity u = 0, 

Dl/Db: <p(u^0) ~ (i + pu 2 + ••• , ip{u -> 0) ~ V (0) + ^ {1) u 2 + • • • , 

D2/DA: 0(u -> 0) ~ /z + pit + ■ ■ ■ , ip(u -> 0) ~ V (0) + H , 

D3/D3: <p(u^0) ~ p + /ilogu + • • • , ^(u -»• 0) ~ + V> (0) log u + ■ ■ ■ . (16) 

Notice that the asymptotic behaviors for (f> and tp between the conformal and nonconformal mod- 
els are different — there is a logarithmic term for the D3/D3 model, which also happens in later 
fluctuation analysis when investigating the electromagnetic response. This difference makes the 
identification of the source and the operator very different from conventional cases. Here, for the 
D3/D3 model, the constant term ip^ is identified as the operator expectation value because it is 
now normalizable with respect to the logarithmal term while ijj^ as the source. Actually, due to 
this fact, the formula for the AC conductivity of the D3/D3 model is also different from the other 
two models. As mentioned before, to make the phase transition a spontaneous symmetry breaking, 
we should impose 

^(°)=0, V {1) ~(0>- (17) 
For the <f> field, the chemical potential fi should be nonzero and the charge density p is a function 



of it. Actually, D3/D3 model has been exhaustively analyzed in the first reference of [14] from the 
field theory point of view. The action for this system is most easily and elegantly constructed in 
(2,2) superspace. Different modes on probe D3-brane have also been studied there. However, one 
main difference from our work is that it also takes into account the £ (internal space in eq. @) 
dependence of flavor U(l) gauge fieldF]. Therefore, the asymptotic behavior for gauge modes near 
conformal boundary is different from our results. 

Near the horizon u = 1, one must have 0(1) = for its norm to be finite while ip should be 
finite there. We then have following Frobenius expansions for <p and ip near the horizon, 

<p(u ~ 1) = ai{u - 1) + h(u - l) 2 + ci(« - l) 3 + di{u - l) 4 + • • • , 
Tjj{ u ~ 1) = a 2 + b 2 (u - 1) + c 2 {u - l) 2 + d 2 {u - l) 3 + e 2 {u - l) 4 + • • • (18) 

where the coefficients bi, etc. can be uniquely determined in terms of a\ and a 2 once plugging 



these expansion into eqs. (13,14,15) and counting by order of (u — 1). In carrying out numerical 



Since we are interested in the lowest energy state, it is reasonable to ignore the £ dependence for flavor gauge field. 
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computations, one first guesses some values for a\ and a<i- Then one uses this near horizon expansion 
to apply a finite-element differential equation solving method. At the boundary u = 0, the condition 
ip(°) = will filter out the wanted solutions. As in M, we will restrict to the solutions where tjj has no 
nodes because on general grounds the solutions with nodes are expected to be thermodynamically 
unstable and have higher energy. 

We plot in FIG. |l| the condensate (O) versus dimensionless temperature T/T c . Explicitly, they 
approach fixed constants as T goes to zero, as is expected for a superconductor. However, the 
expectation values of the condensates near zero temperature are much larger than predictions from 
BCS theory, which should be explained as the strongly coupled feature of holographic supercon- 
ductors. Moreover, our results are similar to higher dimensional counterparts This may be 
understood as one of the universal properties of holographic methods applied to strongly interacting 
superconductors. 
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FIG. 1. The condensate of the 2D p-wave superconductors for the operator O corresponding to D1/D5 
(top-left), D2/D4 (top-right) and D3/D3 (bottom). 

In the mean field theory for the superconductors, the order parameters have them square root 
behavior, 

(O) ~ (T-T c ) 1/2 when T -»• T c . (19) 
12 



By fitting the curves in FIG. ||, we also find such mean field behavior for the condensates in our 
models: 



D1/D5: (O) w (17.8815T C ) 2 (1-T/T C ) 1/2 as T -> T c (20) 
where the critical temperature when expressed in terms of the charge density is T c = 0.0654409p; 

D2/DA: (O) « 15.6435T C (1 - T/T c ) 1/2 when T -> T c (21) 
with T c = 0.0692836p; 

D3/D3: (O) ss 15.8992T C (1 - T/T c ) 1/2 as T -> T c (22) 

where the critical temperature T c is expressed in terms of the chemical potential as T c = 0.0827962/i. 

From these numerical results, we conclude that all three models nearly give the same results 
when concerned with superconducting condensates. 

Before closing this subsection, we plot the results for the free energies for the three systems, 
which can be taken as one evidence for that the superconducting phase transition does happen 



when decreasing the temperature to the critical value T c . With the equations of motion (|13 ,14,1-5) 
for the backgrounds ip and (j), we can reduce the action in eq. (||) to some simpler expression. Then, 
the free energy difference between the normal and superconducting phases is 

AEjv-sc = V f duV^Gg xx g tt [<p(u) V>(u)] 2 = VA_F, (23) 
Jo 

where V is a model dependent factor and we will not give its explicit expression here since it does 
not affect later arguments. FIG. [2] is for the plot of the free energy difference AF. We clearly see 
from these plots that the superconducting phases are thermodynamically favored below the critical 
temperature T c . 



III. ELECTROMAGNETIC FLUCTUATION AND AC CONDUCTIVITY 
A. Fluctuation analysis: formula for the conductivity 

In this section we move on to the study of the electromagnetic response for the superconductor 
models constructed in previous section. In particular, we are concerned with the AC conductivity 
ct{uj) and assume that the fluctuations have no spatial dependence. We focus on the following 
decoupled sector {a^ (t, u), a 2 (i, u), a%.(t, u)} for the gauge field fluctuations and take the Fourier 
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FIG. 2. The dimensionlcss free energy difference between the normal and the superconducting phases for 
D1/D5 (top- left), D2/D4 (top-right) and D3/D3 (bottom), respectively. 



ansatz for them, say {a} (t, u), af(i, u), a%(t, «)} ~ e lujt {aj (ti), af (w), a^(zi)} . Then the linearized 
version of eq. (0) when considering the fluctuations are as follows, 



D1/D5 



i // i i / 

a t a t + 

u 

2" 1 2' 
a* a t 



u 



/(«: 



< + 



/(«) 
1 

u 



(iiutfjal + ^ 2 a 2 ) = 



m 

/(«) 



/ 2 H 



(24) 



where dimensionless frequency For later convenience, we in the below list corresponding 

equations for D2/D4 and D3/D3 models, 



D2/D4 : 



a! + 



.2" 



/(«: 

7R 



0^ = 



(25) 



/(«) * / 2 M 
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D3/D3 



a 



a 



+ 

2" 1 



1 V 1 

~ a t + ~F7 — -i 

U f(u) 



Hal 



9' 

+ -a* 



/(«: 



(zw^ + i> a t ) = 



(26) 



a// 

a' + 



/ 2 H 



[— Cb 2 a? x + ^t/jaj + iutpaf] = 0, 



where w = JfW- Actually, these modes are not independent because the radial gauge aj = 
has been chosen in deriving these fluctuation equations. This gauge choice gives two constraint 
equations for this decoupled sector. We do not present them here as these constraints have no 
effect on the definition as well as the numerical computations for the conductivity, which will be 
clear later. The AC conductivity is defined by following Kubo's formula, 



ICO 



(27) 



Therefore, our aim to produce the electromagnetic response of the superconducting models we 
constructed in section [ij is then reduced to calculation of the retarded Green's function Gr(uj). 
Under gauge/gravity duality approach, a good prescription for the retarded Green's function can 



be found in [31]. 

For the retarded Green's function, we should impose ingoing wave boundary condition at the 
horizon for these fluctuation modes. A simple Frobenius analysis for eqs. ( p4| , ^5|j2^ ) near the 
horizon reveals that 





-u) a 




(l -„) + (£ w (l - 


u? + af\l- 


-uf + ---~ 




< 4 = (i 


-u) a 




u) + a\ (2) {1 - uf 


+ af\l-uf 


+ ••• 




(28) 


ai = (1 


-u) a 






+ af\l-uf 


+ ••• 







where we have used the linearity of eqs. to set the scale of a x to 1 at the horizon. The 

indices appearing in these equations are —iui/6, —ioj/b and —iCo/A for the three models, respectively. 
The coefficients in above equations can be uniquely determined once plugging these expansions into 



eqs. (^^ |_J) an d counting by powers of (1 — u). Therefore, eqs. (p§|) can provide initial conditions 
for these second order differential equations. We in fact use these power solutions to do numerical 
integration from the horizon to the conformal boundary by the mathematica NDSolve. 

At the AdS conformal boundary u = 0, the general solution to the equations of motion is 
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distinguished between the D1/D5 (D2/D4) and D3/D3 models. More specifically, 



D1/D5 (D2/D4) : 



,1(0) 



and 



D3/D3 



.4 



.4 



Ar+Ar>u\ U ) + 

Af ] +A 2 t {1) u 2 (u) + 
Af ) +Af ) u\u) + 



)) +4 (1) logn + 
2(0) + A t 2(1) logn + 

,3(1) 



(29) 



(30) 



4 = Al {0) +A%^> logu + 



where we have represented the boundary expansions for the D1/D5 and D2/D4 models together 
and the expressions in the parenthesis are for the D2/D4 system. Explicitly, the logarithmic terms 
appear again in the boundary behavior of the gauge field fluctuations for the D3/D3 model. 

As argued in ||], the conductivity is a physical quantity and should be gauge invariant. We 
should construct a new mode from a\,a 2 , a x and this mode should be invariant under the gauge 
transformation that keeps our gauge choice. The details for the constructions of the gauge invariant 
modes can be found in the original paper and we in the following write down this mode directly, 



iua 2 + 4>a\ 
6 2 — Cj 2 



(31) 



Plugging the boundary behavior in eq. (29) into the newly defined mode in eq. (|3l| ) and expanding 
it near u = results in 



D1/D5(D2/D4) : of = A 3 X (0) + V(u) + 



(32) 



with 



<°) = 4SW> jpV> = 43(1)^,(1) 



-/1 5. — A 



i) , „,.(D ^A 2(0) +^A\ 



(0) 



A 6K ' -i- 7/4 

x i J *x ^x ' T ,.2 ~2 



(33) 



Then, the formula for the conductivity can be straightforwardly written down for D1/D5 and 
D2/D4 systems, 



±4? 



(34) 



Notice that the above formula is the same as the equation (4.19) of We can expect that 
numerical results for the conductivity of these two models should have some similarities with the 
results reported in H for a xx and this does happen for our results. 



16 



Involving the D3/D3 model, we have the corresponding results for the mode a| near u = 0, 

D3/D3 : a% = A 3 X (1) + log u + • • • , (35) 

with 

i3 W _ A 3 (0) i3 W _ ^3 (1) ,3^ 
^X — ^x ) ^x — ^x t ■ \ ou ) 

Due to the appearance of the logarithmic term log u in the asymptotic behavior for the fluctuation 
modes, we should identify A% as the source and A%} ^ as the expectation value of the dual 
operator. Moreover, the definition for the retarded Green's function as well as the conductivity 
should be modified to 



i3 (1) 1 A* {1) 

G H ( W ) = __2- and a( W ) = -^. (37) 



We have seen that this formula is greatly different from the counterparts eqs. (33,34) for D1/D5 
and D2/D4 models. In the next subsection, we will see that this difference will be reflected in the 
imaginary parts of the conductivity. 

B. Numerical results for the conductivity 

We report our numeric results for the AC conductivity for all three models in this subsection. 
Before this, we have a brief explanation for our numerical method. We use the horizon expansions 
as in eq. ( |28| ) to generate initial conditions for the equations of motion for the fluctuations. We 
then use the mathematica NDSolve to numerically solve these eqs. (|24|, 25 , 26|) . The expansion 



coefficients appearing in the conductivity formula can be directly read off from boundary behavior 
of different modes as in eqs. (p9|,30) once the backgrounds ip and <p as well as the frequency to are 



specified. 

In FIG. (QODflO D'- wc P lot the real and 

imaginary parts of the AC conductivities for all the 

three models. 

We conclude with some explanations of these figures for the conductivity: 
Firstly consider the normal phase, i.e. T/T c = 1. For the D1/D5 model, Re<r(a;) is a constant 1 
and Imff(w) is always zero. This is consistent with higher dimensional models. On the other hand, 
for D2/D4 and D3/D3 models, Recr(u;) decreases monotonously as increasing the frequency uj/T. 
Moreover, Hea(co) of D3/D3 approaches to zero at high frequency while Recr(w) for D2/D4 goes 
to a nonzero constant. This feature has also been found by studying the flavor U(l) dynamics of 



D3/D3 model or BTZ black hole in ||, |TJ, |7|. 
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FIG. 3. The real part of the AC conductivity with different condensates corresponding to T/T c = 
1.0,0.603055,0.277291,0.193113,0.126896 for D1/D5 from top to down. 



FIG. 4. The real part of the AC conductivity with different condensates corresponding to T/T c = 
1.0,0.464611,0.254859,0.141036,0.0978091 for D2/D4 model from top to down. 

When decreasing the temperature, we go to the condensed phase, i.e. (O) ^ 0. We find that 
Rea(oj) for all three models display gap formation, which is similar to the findings in [||, || [l7| . 
More specifically, Re<r(u;) is very small in the infrared (i.e. for small uj/T in the figures) and 
then rises quickly at some critical frequency ui g . Near u = 0, there appears a delta peak in the 
conductivity, which can be explained by the Kramers-Kronig relations. This relation relates the 
real and imaginary parts of casual quantities and can detect the distributional parts of them. Take 
the D1/D5 system as an example. The Kramers-Kronig relation states that, 



where P denotes the principal part of the integration. It is clear from this formula that Re(u;) 
has a delta peak, Recr(w) ~ only when Ihkt(w) has a pole at u = 0, Re<7(w) ~ 1/w, and 

vice versa. The delta peak at lo = indicates a DC superconductivity. Similarly, the pole in the 
imaginary part of the AC conductivity at finite oj can also be understood from the Kramers-Kronig 




T 




(38) 
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FIG. 5. The real part of the AC conductivity with different condensates corresponding to T/T c = 
1.0(Solid),0.664919(Dashed),0.269669(Dotted),0.I89074(Thick Dashed) for D3/D3 model. 
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FIG. 6. The imaginary part of the AC conductivity with different condensates corresponding to T/T c = 
1.0,0.603055,0.277291,0.193113,0.126896 for D1/D5 from left to right and top to down. 

relation. As mentioned before, this pole can be easily seen from the formula for the conductivity 
eq. (|33|). From eq. (|38|), we can conclude that a simple pole in Im[cr(w)] at lo = loq implies a 
delta-function contribution 5(cj — uiq) to Re[a(u>)]. We already see that this peak becomes more 
higher and narrower as decreasing the temperature. However, as for as we understand, we do not 
think of it as a massive excitation. We rather take it as artificial as we use the newly defined mode 
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FIG. 7. The imaginary part of the AC conductivity with different condensates corresponding to T/T c = 
1.0,0.464611,0.254859,0.141036 for D2/D4 from left to right and top to down. 

to plot the conductivity, which in some sense implies that it is not a good choice of gauge invariant 
mode used in this work and the pole is introduced by hand, however inevitable as far as we know. 

In the superconducting phase, there is a pole in Im<r(u;) for D1/D5 and D2/D4 models, which 
is in contrast to the D3/D3 model. Actually, from the formulae of the conductivity for all three 
models, we have seen some differences between them and the peak should appear at ujq = f/, from 
eq. (|3^) with the residue proportional to the condensate (O). However, for D3/D3 model, we can 



easily read off the uj g if we follow the analysis of [32] and define oj 9 as the frequency which minimizes 
Im(o;), 

T 

-± « 25 when — = 0.269669. (39) 

When taking large oj limit, all the results go to those of the normal phase, which is expected 
from general grounds, the large u will wash out the effect of the superconducting condensates. 



This can also be seen clearly from the eqs. (p4| : 25,2 



With these comments in mind, we can conclude that all three models can reproduce some basic 
features of the 1+1 dimensional p-wave superconductor, like the DC infinite conductivity (a signal 
for superconductivity) and gap formation when decreasing the temperature. The results of D1/D5 
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FIG. 8. The imaginary part of the AC conductivity with different condensates corresponding to T/T c = 
1.0,0.664919,0.269669,0.189074 for D3/D3 from left to right and top to down. 

model are more related to the higher dimensional counterparts. The D3/D3 model for p-wave 
superconductor have some common feature with its s-wave case as studied in g7[- 

IV. SUMMARY 

In this paper we have taken the D-brane probe approach to explore some properties of holo- 
graphic p-wave superconductor in 1+1 dimensional spacetime. In the large N limit, we can bypass 
the no-go theorem (the CMW theorem) which forbids the superconducting phase transition. We 
found that all the three models are quantitatively similar in producing some key features of the 
p-wave superconductor, like the mean field behavior of the superconducting condensate near T c , 
the DC delta peak and the gap formation, etc. These are common with the higher dimensional 
superconductors under the holographic approach. Therefore, these characteristics can be regarded 
as universality of the holographic method. Besides these, we also found some particular ones for 
different models, especially for the AC conductivity. Specially, the high frequency limits of the con- 
ductivities are different: for D1/D5, it goes like the higher dimensional counterparts, Re<7(w) — > 1 
and Im<r(u;) — > 0; for D2/D4, both Recr(u;) and Imcr(u;) approaches some nonzero constants but 
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very slowly compared to the other two models; the conductivity of the D3/D3 model is basically 
the same as its s-wave counterpart, which is first studied in fllTf] . However, the condensates found 
here is greatly different from the results of [17] where (O) decreases or increases in the zero tem- 
perature limit. Leaving aside these model dependent differences, we can conclude the holographic 
approach to the superconducting phase transition can reproduce some common features and give 
us some directions for strongly correlated condensed matter system. 

Actually, we in this work only explored some basic aspects of the 1+1-dimensional supercon- 
ductor from holographic viewpoint. There are some interesting problems deserved further investi- 
gations. The first one is to go beyond the zero quark mass limit and study its effect on the super- 
conducting phase transition and conductivity. There is only one conductivity for 1+1-dimensional 
spacetime and it looks like a xx (uj) in the higher dimensional case. However, such quantity in color- 
ful superconductors from probe D-brane approach has not been computed. So we hope to come to 
this question in the near future and have a complete comparison between our work and its higher 
dimensional counterparts. 
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